A family of sure-success quantum algorithms for solving a generalized Grover search 

problem 
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This work considers a generalization ol Grover's search problem, viz., to find any one element in a 
set of acceptable choices which constitute a fraction / of the total number of choices in an unsorted 
data base. An infinite family of sure-success quantum algorithms are introduced here to solve this 
problem, each member for a different range of /. The nth member of this family involves n queries 
of the data base, and so the lowest few members of this family should be very convenient algorithms 
within their ranges of validity. The even member A^n of the family covers ever larger range of / for 
larger n, which is expected to become the full range < / < 1 in the limit n — > oo. 

PACS numbers: 



Quantum computing as a new powerful approach to 
solve difficult computational problems is still in its in- 
fancy. Only a handful useful algorithms have been pro- 
posed so far. Most of them fall into two categories: Those 
for factorizing large integers, and those for "searching a 
needle in a haystack" , or finding the only acceptable el- 
ement in a large unsorted data base. The main idea in 
the former category is due to Shor, [|| and in the lat- 
ter category is due to Grover. ||] Here we wish to con- 
sider a generalization of Grover's search problem, viz., 
to find any one element in a set of acceptable choices 
which form a fraction / of the total number of choices 
in an unsorted data base of size N. ^ An infinite fam- 
ily {An \ n — 1,2,3,- ■ ■} of qunatum algorithms is intro- 
duced here, each similar to one stage of Grover's algo- 
rithm, except that, unlike Grover's original algorithm, 
which requires iteration to some optimal stage, which in 
general is still not a sure-success algorithm, [Q here each 
member is an independent sure-success algorithm within 
its range of validity. Each member of the family intro- 
duced here is characterized by an iteration number, n, 
in the sense introduced in the original Grover algorithm. 
This number is also the number of times the data base is 
queried. Here we only analyze four members of this fam- 
ily, corresponding to the iteration numbers 1,2, 4, and 6. 
We find that Ai is vahd for 0.25 < / < 1.0; A2 is valid 
for 0.095491502 •••</< 0.65450849 ■ ■ ■; A4 is valid for 
0.030153689 •••</< 0.88302222 • • •; and Ae is valid 
for 0.014529091 •••</< 0.94272801 • • •. These results 
strongly indicate that by using A2n of ever larger n, an 
ever larger range of / can be covered which in the limit 
of n ^ 00 approaches the full range < / < 1, but a 
general proof has not yet been obtained. The validity or 
non- validity of this statement, and the properties of the 
odd members of this family, will be discussed in a future 
work. All members of this family of algorithms are char- 
acterized by two phase parameters, 9 and </>. These two 
parameters are individually adjusted in order to make 
each member a sure-success algorithm. I find that (i) 



at least for the members A2, A4, Ae, but most-likely 
also for all higher even members of the family, (j) = 29 
is an acceptable choice for (j), (ii) for each of the even 
members A2, A4, Ae, and most-likely also for each of all 
higher even members of the family, the required value for 
6 for it to work is a unique function of / just inside the 
boundary of its validity range of /, but the number of 
acceptable values of 9 graduately increases to n deep in- 
side the validity /-range for A2n ■ The algorithm member 
Ai, on the other hand, requires (j) = —29, then 9 depends 
uniquely on / within the validity range. No other odd 
members have yet been analyzed. In all cases studied, I 
find the required 9 and (j) to be independent of N, and 
to only depend on /. There is strong indication that this 
statement is true for all members of the family. 

All members of this family {An | f^ = 1, 2, 3, • • •} are 
achieved with two unitary operators which generalize the 
two corresponding operators introduced by Grover: 

In a Hilbert space spanned by a set of N orthonormal 
basis states {\i > | n = 1, 2, 3, • • • , N}, each of which rep- 
resents one element of the data base, Grover introduced 
an unitary operator, which I shall denote as F^, which 
changes the sign of the i^'s amplitude Ci, in any quan- 
tum state l^ii >= X^i^i I* >■ This operator is gener- 
alized to the operator F^"'' , which introduces the extra 
phase factor — e*"^ to each of the amplitudes {Ci, | G a}, 
where a denotes the set of acceptable elements in the 
data base. j|] Mathematically, 
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where / = X^il^i K *l is the identity operator. For 
(j> — 0, and a containing only one element v, this operator 
reduces to the operator F^, introduced by Grover. 

A second unitary operator introduced by Grover is the 
"inversion about the mean" operator, which can be writ- 
ten in the form: 
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j(2/iV)-4,]|z><j|. 



(2) 
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I generalize it to 

Oe = ^[(2cos0/7V) - e^'S.^j] \i >< j\ , (3) 

which reduces to Grover's "inversion about the mean" 
operator ii 9 ~ 0. That Oe is unitary can be easily 
verified. It is also easy to show that it is the most general 
unitary operator of the form J2i jii^ + I* 
if one disregards an unimportant overall phase factor. I 
am not aware of any earlier published work introducing 
this umitary operator. 

Since F^"^ and Og are both complex operators, I also 
need their hermitian conjugate operators, i^j""*^ and Og, 
which are also the inverse operators of F^"-* and Og, re- 
spectively. Actually they are simply F^J and O^g. 

Before any algorithm is applied, every element in the 
data base should be regarded as to have equal probability 
of being the right choice. Grover represented this fact by 
starting with the quantum state: 

N 

1*0 >= (i/Vlv)^K>, 

i 

i.e., the state with every Ci — so that the 

probability of finding any element of the data base is 
|Cip = The quantum algorithm he introduced is 

to repreatedly apply the unitary operator product OFi, n 
times on the state |5'o >, followed by a measurement to 
cause the state to collapse to one of the basis states. He 
showed that when n is of an optimal value of the order of 
•\/iV, All |Cip will be very close to zero except the par- 
ticular one |Cyp, corresponding to the desirable element 
v in Grover's search problem, which will be very close to 
unity. However, except for some special values of TV, one 
will not obtain exact unity for |Cj/p, and exact zero for 
all other |Cip. Thus Grover algorithm is in general not 
a sure-success alorithm, even in theory, when potenial 
implementation errors are not taken into account. We 
generalize Grover's algorithm to a family of sure-success 
algorithms, each member of which is characterized by an 
integer n. Denoting these member algorithms as {^n}, 
then the even [(2n)th] member {^2n} are defined as ap- 
plying the unitary operator product A = OJf^"^ ^OgF^"'^ 
n times to the state 1^*0 >, followed by the same measure- 
ment used in the Grover algorithm. The odd [(2n-|- \)th] 
member {-42n+i}, is to apply the unitary operator prod- 
uct OgP'^^'^k^ to the state >, before the same mea- 
surement is made. Thus {.4„} makes n queries of the 
data base. "Sure success" of each of these algorithms is 
achieved by adjusting the two parameters 6 and cj) so that 
all |Cip, with i not belonging to the set a of the general- 
ized Grover search problem introduced here, are exactly 
zero. All |Cip with z e a will then be exactly equal to 
1 / {fN) , where fN = Na is the number of elements in 



the set a, since probability is conserved by unitary oper- 
ations. Below we show how this is done explicitly for the 
four members Ai, A2, A4, and Aq. After that I will spec- 
ulate about all even members A2n of the family, leaving 
the odd members higher than the first to be discussed in 
a later work. 

Consider first the algorithm member Ai . One has the 
identity: 

OeF^^^l^o >= [2 cos 0(1 - / - /e^^) - e'^F^'^'H^o > . 

(4) 

Since the operator is equivalent to an identity op- 
erator in the subspace corresponding to all unacceptable 
elements of the data base, sure success of this algorithm 
is achieved by demanding 

2cos6l(l-/-/e"^)-e*'' = 0. (5) 

which has the solution <j) ~ —26', and 

0=(l/2)cos-i[(l/2/)-l]. (6) 

Note that if {(p, 9) is a solution, then {—(j>, —9) is also 
a solution. This is true for all higher members of the 
family also, and one can easily see why. Equation (|^) 
has solution only for 1/4 < / < 1, which is the va- 
lidity range of this algorithm. Within this range, I 
have plotted 6* as a function of / in Fig. |l| assuming 
9 > Q. The following special cases are of interest: (i) 



1.2 




0.2 0.4 0.6 0.8 1 

f 

FIG. 1. Plotted is 9 versus / for the algorithm ^1. 

For / = 1/4, I find (p — 9 — and the operators re- 
duce to those introduced by Grover, and this algorithm 
becomes a special case of Grover's algorithm, (ii) For 
/ = 1/3, I find </) = ±7r/3 and 9 = T7r/6. (iii) For 
/ = 1/2, I find (/) = ±7r/2, and 9 = T7r/4. (iv) For 
/ = 2/3, I find <?!) ±104.477- • = ±0.580430.. .tt and 
9 = =f52.2387---° ^ =f0.290215 • • -vr. Finally, (v) for 
/ = 1, I find (f) = ±27r/3, and 9 = =F^/3, but in this case 
the operator product OeF^"' acting on |$o > simply re- 
produces |<i>o >■ 
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Next, let us consider the second member A2- One has 
the identity: 



A|<i>o > =mcosef\{l-f-fe 



-(2cos0)e-^«(l-/-/e^^)Fi")^}|<i>o> 
(i)t^ 



(Ai-SiFr')|<f>o> 



(7) 



(Note that Ai = \Bi\ 



) Thus to ensure that this 



is a sure-success algorithm, one needs only demand Ai — 
Bi = 0. The imaginary part of this condition can be 
written as 



Im( Ai - Si ) (2/ cos 9) [sin(0 -9)- sin ( 



0. (8) 



so it can be satisfied with </> = 29. (It is easy to see that 
cos 9 7^ 0.) Then the real part of this condition reduces 
to 



Re(Ai - Si) = 1 + Af^? - 16/(1 - f)^i^ = 
where jj. = cos 9. It has the solution 



^cos-{^(^[y|-3+(4/-3)]}. 



(9) 



(10) 



This equation has solution only if 0.095491502 • • • < 
/ < 0.65450849 • • •. Within this range, I have plotted 9 
as a function of / for this algorithm in Fig. ||, assuming 
6* > 0. 




FIG. 2. Plotted is 6 versus / for the algorithm A2- 

Note that with the lagorithm A2 we can cover / down 
to slightly below 0.1. 

Next, let us consider the algorithm member ^4, leav- 
ing ^3 and higher odd members for future discussion, 
since they are deemed less important. I have first es- 
tablished the following theorem: If A^j^o >= [^n — 
B„Fi'^)t]|«'o >, then 
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l-o > = {[AiA„ 



-2ie n* 



-[B,A^-e-^''Br,]F^''^^\^o> ■ (H) 



That is, 



An+l 
Bn+1 



Ai -B*e-2'6 
Bi 



An 
Bn 



(12) 



Thus A2 = - [2cos(26') + e2*«]|Bi|2 + e'^'" and 

B2 — [\Bi\^ — 2 cos(20)]_Bi. To ensure sure-success for 
this algorithm, one needs to require A2 — B2 = 0. It is 
easy to show that 

Im(A2 - S2) = [|Si p - 2 cos(26')]Im(Ai - Bi) . (13) 

I shall consider in a future work the possibility of satisfy- 
ing this equation by setting the first factor equal to zero. 
Here I concentrate on the fact that due to its second fac- 
tor this equation can be satisfied by letting (j> = 29. Then 
9 is given by 



Re(^2 - B2) = 1 + 8/^2 _ 48^(1 _ j)^4 
-64/2(1 - /)// + 256/2(1 - /)2/ = . 



(14) 



I have plotted as a function of / for this algorithm in 
Fig. H assuming 9 > 0. It is seen that solution exists only 




FIG. 3. Plotted is 6 versus / for the algorithm Aa- 

for 0.030153689 •••</< 0.88302222 • • •, and that in 
the narrower range 0.25 < / < 0.58682408 • • • A second 
solution for 9 appears for each /. It should be obvious 
that this algorithm is valid for those values of / only, for 
which at least one solution for 9 exists, thus the larger / 
range is also the validity range of this algorithm. 

Finally, let us consider the algorithm member Ae- 
Eq. (H) allows me to obtain A3 = \Bi\^ - [4cos(26') -I- 



e^'"]\Bi\^ + 2[cos(46l) + 1 



and B3 



- 4cos(26l)|Bi|2 + [2cos(46l) + l]}Bi. Thus I find 



lui{As - B3 



{[\Bi\' 



2cos(26i)]^ 



l}lm{A,-B,). 

(15) 



Again, I shall not consider here letting the first factor 
equal to zero. Then again (f) = 29 from Im(yl3 — B3) = 0, 
and 9 is given by 
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Re{As - Bg) = 1 + 12/^^2 _ 96/(1 - /)/ 
-256/2(1 - /)/i6 + 1280/2(1 - fff/ 
+ 1024/3(1 - f)^fi^° - 4096/3(1 - /)3//2 = . (16) 

I have plotted as a function of / for this algorithm in 
Fig. assuming 9 > 0. It is seen that solution exists only 



is then the most convenient, since it requires the least 
number of queries of the data base. 
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FIG. 4. Plotted is 9 versus / for the algorithm Ae- 

for 0.014529091 •••</< 0.94272801 • • •, which is the 
validity range of this algorithm. In the narrower range 
0.12574462- •• < / < 0.78403237 •• • a second solution 
for 9 appears for each /, and in the even narrower range 
0.32269755 •••</< 0.56026834 • • • a third solution for 
9 appears for each /. 

A trend is clearly established by the above study of 
the first three even members. It strongly suggests that 
for all even members, (i) (j) — 29 is always a valid solution, 
with 9 depending on /, but not on N; (ii) the /-range 
in which at least one 9 value exists becomes ever larger 
if A2n of ever larger n is considered, with the n — > oo 
limit being very likely the full range < / < 1; (hi) in 
general the number of valid choices for 9 increases to n 
deep inside the validity /-range for A2n- General proofs 
of these statements have not yet been obtained. 

In summary, an infinite family of sure-success quantum 
algorithms is introduced here for solving the generalized 
Grover search problem of finding any one element of a set 
of acceptable choices which constitute a fraction / of all 
elements in an unsorted data base. This is achieved by 
two unitary operators each containing a phase parame- 
ter. These operators are generalizations of the two opera- 
tors introduced by Grover for his original search problem. 
The two phase parameters are adjusted for each member 
of the family to ensure its sure-success, which is found 
possible only within a different /-range for each mem- 
ber of the algorithm family. An infinite sub-family (the 
"even" members) appears to have the property that the 
validity /-range of a lower member is totally embedded 
inside that of a higher member, with the limit being very 
likely the full range < / < 1. As long as / is within 
the validity range, the lowest member of the sub-family 
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